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We study CP violation and the contribution of the strong pion-pion 
interactions in the three-body B^ — > tt^tt^tt^- decays within a quasi two- 
body QCD factorization approach. The short distance interaction ampli- 
tude is calculated in the next-to-leading order in the strong coupling con- 
stant with vertex and penguin corrections. The meson-meson final state 
interactions are described by pion non-strange scalar and vector form fac- 
tors for the S and P waves and by a relativistic Breit-Wigner formula for 
the D wave. The pion scalar form factor is calculated from a unitary rel- 
ativistic coupled-channel model including tttt, KK and effective (2tt)(2tt) 
interactions. The pion vector form factor results from a Belle Collabora- 
tion analysis of r~ —y tt~tt q v t data. The recent B^- —y 7r ± 7r T 7r ± BABAR 
Collaboration data are htted with our model using only three parameters 
for the S wave, one for the P wave and one for the D wave. We find not 
only a sizable contribution of the S wave just above the tttt threshold but 
also under the p(770) peak a significant interference, mainly between the 
S and P waves. For the B to /2(1270) transition form factor, we predict 
F B f 2 (m^) = 0.098 ± 0.007. Our model yields a unified unitary descrip- 
tion of the contribution of the three scalar resonances /o(600), /o(980) and 
/o(1400) in terms of the pion non-strange scalar form factor. 
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PACS numbers: 13.25.Hw, 13.75Lb 

1. Introduction 

Three-body charmless hadronic B meson decays offer one of the best 
tools for studies of direct CP violation and provide an interesting testing 
ground for strong interaction dynamical models. The present work, part of 
a program devoted to the understanding of rare three-body B decays [U [21 
El 13] , is motivated by the recent BABAR Dalitz-plot analysis of the B^ — > 
tt^tt^tt^ 1 decays [S]. In an isobar model description, the authors of Ref. [5] 
find evidence for the /o(1370) but, within the current experimental accuracy, 
no significant signal for the /o(980). The /o(600), not explicitly included in 
that analysis, could be part of the non-resonant background. Furthermore, 
there is a small but visible contribution of the /2(1270) resonance [5]. 

Here, the aim is to provide a phenomenological analysis of the S ± — > 
tt^tt^tt^ decay channels relying on the QCD factorization scheme (QCDF) 
in the tttt effective mass range from threshold to 1.64 GeV. The focus will 
be set on the final state tttt interactions involved since a partial wave anal- 
ysis of the Dalitz plot should use theoretically and phenomenologically well 
constrained tttt amplitudes. 

Studies of B decays into two-body and quasi-two-body final states have 
been performed in the QCDF framework [61 El O \TU[ \TT\ I12j . The naive 
factorization approach is a useful first order approximation which receives 
corrections proportional to the strong coupling constant a s at scales m& 
and y / A-QCDmb and in inverse powers of the b quark mass ra\, [13]. In the 
present study, we propose an extension of these results to the three-body 
decays B^ — > 7r ± tt + tt^. 

The role of the /o(600) (or a) in charmless three-body decays of B 
mesons has been examined by Gardner and Meifiner [8] in B° — > tt + tt~tt° 
decays. Within QCD quasi two-body factorization approach their /o(600)7r 
amplitude is described by a unitary pion scalar form factor constrained by 
tttt scattering and chiral dynamics. This is different from the relativistic 
Breit-Wigner parametrization used in most experimental analyses and in 
some theoretical studies, for example in |14j . This has led to improved 
theoretical predictions; the contribution of the /o(600)7r channel has been 
found to be important in the range of the dominant 

A 

intermediate state. 

However, in recent B° — > tt + tt~tt Dalitz plot analyses |151 1 16] no contribu- 
tion from B° — > /o(600)-7T° channel has been found. This could be linked to 
the present limited statistics in the low effective tttt mass region. Further- 
more, such a contribution could be hidden in the nonresonant amplitude 
introduced in the experimental analysis. Nevertheless we will show that the 
contribution of the S wave is important in the B ± — > 7r ± tt + tt decays. 
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Charmless three-body decays of B mesons have also been investigated 
by Cheng, Chua and Soni [12] in the framework of quasi two-body factoriza- 
tion approach using resonant and non-resonant contributions. In particular 
they have calculated the B~ — > ir + ir~Tr~ branching fractions and CP asym- 
metries and found a small rate for B~ — > /o(980)-7r~decay. 

An achievement in the theory of B decays into two mesons is the confir- 
mation of the validity of factorization as a leading order approximation. No 
proof of factorization has yet been given for the B decays into three mesons. 
However, three-body interactions are suppressed when specific kinematical 
configurations with the three mesons quasi aligned in the rest frame of the 
B meson are considered. This is the case in the effective 7r + 7r~ mass region 
smaller than 1.64 GeV in the Dalitz plot where most of the tt + tt~ resonant 
states are visible. Such processes will be denoted as B^ — > tt^I^tt'}, the 
mesons of the [7r + 7r~] pair moving more or less, in the same direction in the 
B rest frame. Then, it seems reasonable to postulate the validity of factor- 
ization for this quasi two-body B decay [T7] assuming that the [7r + 7r~] pair 
originates from a quark-antiquark state. 

In the factorization approach the B^ — > ir^ [vr^7r^] decay amplitudes 
are expressed as a superposition of appropriate effective QCD coefficients 
and two products of two transition matrix elements. The transition ma- 
trix elements between the B^ meson and the 7r^ pion multiplied by the 
transition matrix elements between the vacuum and the [^7^] pion pair 
correspond to the first of these products. Here, in the tt^ tt^ center of mass 
frame, the bilinear quark currents involved force the [vr^vr^] pair to be in 
S or in P state. The second term is associated to the product of the tran- 
sition matrix elements between the B^ meson and the [vr^"^] pion pair 
in S, P or D state by the transition matrix elements between the vacuum 
and the pion. The [^2 7T 3]s,P transition matrix elements to the vacuum 
are proportional to the pion scalar and vector form factors. We assume 
that the B^ — > ^2^3 ma t r i x elements are expressed as products of the 
B^ — > [t^2 7: 3]s,p,d transition form factors by the relevant vertex function 
describing the decay of the [^2 7T 3~]s,p,d state into the final pion pair. The 
vertex functions are in turn assumed to be proportional to the pion scalar 
form factor for the S wave, to the vector form factor for the P wave and to 
a relativistic Breit-Wigner formula for the D wave. Here, a single unitary 
function, namely the pion non-strange scalar form factor, describes then the 
three scalar resonances, /o(600), /o(980) and /o(1400) present in the 7r + 7r~ 
interaction. 

In Sec. [2] we present the model used in the analysis. Sec. [3] is devoted to 
the construction of the pion scalar and vector form factors. The pertinent 
observables and the fitting procedure are described in Sec.[3]while the results 
are discussed in Sec. [SJ A summary and some perspectives are outlined in 
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the final Sec.[6j The detailed derivation of the decay amplitudes is presented 
in the Appendix [A] while Appendix [B] gives the system of equations to be 
solved to obtain the parameters fixing the low-energy behavior of the pion 
scalar form factor to be that of one loop calculation in chiral perturbation 
theory. 

2. Decay amplitudes 

The amplitudes for the non-leptonic decays of the B meson are given as 
matrix elements of the effective weak Hamiltonian [61 [7] 

H eff = -| A p [Ci(%+C 2 (%+ , £c i Oi+C 7l 7l +CsgOsg\+h.c., (1) 

* p=u,c i=3 

where 



An = V ub V: d , X c = VM, (2) 

the Vppi (p' = b, d) being Cabibbo-Kobayashi-Maskawa quark-mixing ma- 
trix elements. For the Fermi coupling constant Gj? we take the value 
1.16637 x 1(T 5 GeV" 2 . The d(p) are the Wilson coefficients of the four- 
quark operators Oi(fx) at a renormalization scale [i. The 0\ 2 are left-handed 
current-current operators arising from VF-boson exchange, Oj=3_io are QCD 
and electroweak penguin operators involving a loop with auorc quark and a 
W boson, O-j^ and 0% g are the electromagnetic and chromomagnetic dipole 
operators [7J. 

Let pb be the four- momentum of the meson and p\ that of the 
isolated ir^. Let then p2 denote the four-momentum of the ir + and p% that 
of the 7r~ of the interacting [7r + 7r~] pair in the B rest frame. One has 
Pb = Pi + P2 + P3 and we introduce the invariants Sjj = (pi + pj) 2 for 
i,j = 1,2,3 with i < j. For the B~ — > -ir~ [tt + tt~]s,p,d amplitude, we 
work in the center of mass frame of the tt + tt~ pair of pions with respective 
four-momenta p2 and p% (or p\ and p2 for the symmetrized amplitudes). 
These two pions will be either in a relative S, P or D state. In the following 
we derive the amplitudes for the B~ — > ir~ [tt + 7t~]s,p.d processes. The 
transcription to the B + — > tt + \n + n~\s p d processes is straightforward. 
Applying the QCD factorization formula for the B~ — > ir~ [t^ + ^]s,p,d 
process, the matrix elements of the effective weak Hamiltonian ([T|) can be 
written as [7J 

(vr~Oi) \TT + (p 2 )'K'{pz)\s,P ) D\H e ff\B'{p B )) = 

%Y. X p Ik + *-]s,p,d\T p \B-), (3) 

* p=u,c 
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to which must be added the symmetrized term 

(K~(p3)[TT + (p2)Tr~(pi)]s,p,D\H eff \B-(p B )). With M 1 = tt~ and M 2 = 
[ir + Tr~]s,p or Mi = [it +/ k~}s,p.d while M 2 = tt~ , one has 



(ir~ [ir + -K^} S) p,D\T p \B~) = (ir~ [k +/ k~]s,p,d\ 

I ai{MiM 2 )5 pu (ub) V -A <8> (du) V -A 

+ a 2 (MiM 2 )(5p U ((i6)y_ J 4 ® (uu) v ~A 

+ a 3 (MiM 2 ) ^(d6) y -A ® (qq)v-A 
q 

+ a^(MiM 2 ) ^(g6)y„ A ® (A/)v-a 

+ ag(MiM 2 ) ^(-2)(<?6) sc _ ps 8) (dg) sc+ps 

3 

+ a 7 (Mi¥ 2 ) ^(J6)v-a ® -e g (qq)v+A 
q 

3 

+ a|(MiM 2 ) J2(-2)(qb) sc - ps (g> -e q {dq) sc+ps 
q 

3 

+ a 9 (Mi¥ 2 ) ^(J6)y_A ® -e q (qq) V -A 
q 

+ a? (MiM 2 ) ^(g6) y -A ® -e g (Jg)y_A} 
I A"), (4) 
where a? are effective QCDF coefficients. 

In Eq.flID, (qiq 2 )vTA = ?n)t(l =F 7s)<?2, (<M2)sc±p S = <7i(l ± 75)<?2 and 
e g denotes the electric charge of the quark q in units of the elementary 
charge e. The sum on the index q runs over u and d and the summation 
over the color degree of freedom has been performed. The notations sc and 
ps stand for scalar and pseudoscalar, respectively. 

At next-to- leading order (NLO) in the strong coupling constant a s , the 
general expression of the a? quantities in terms of effective Wilson coefffi- 
cients is [U] 



G 



a: 




+—H l {M 1 M 2 ) +Pf{M 2 ) 



+ 



Ci±l CFOi s 



N C 47T 



[Vi(M 2 ) 



(5) 



where the upper (lower) signs apply when the index % is odd (even), Nc is 
the number of colors, N c = 3 and C F = (JVg - 1)/2N C . The sums over the 
color degree of freedom have been performed in Eq. Q. Note that in the 
leading-order (LO) contribution Ni(M 2 ) = for M 2 = [7r + 7r~]p and i = 6,8, 
otherwise Ni(M 2 ) = 1. The NLO quantities Vi(M 2 ) arise from one loop ver- 
tex corrections, Hi(M\M 2 ) from hard spectator scattering interactions and 
Pf(M 2 ) from penguin contractions. Here the meson M 2 is the meson which 
does not include the spectator quark of the B meson. The superscript p 
in aP i {MiM 2 ) is to be omitted for i = 1, 2, 3, 5, 7 and 9 since the penguin 
corrections are equal to zero in these cases. In our calculation we shall not 
include the NLO hard scattering corrections nor the annihilation contribu- 
tions which require the introduction of four phenomenological parameters to 
regularize end point divergences related to asymptotic wave functions |9]. 
Although we are aware that such contributions might be important, this 
would bring, at this stage of analysis, too many free parameters. 

In Eq. ([3D the symbol <%> indicates that the different components of the 
matrix elements (tt~ [tt + h~]s,p ) d\T p \B~) are to be calculated in the factor- 
ized form, 



(vr (pi)[7r + (p 2 )vr (jpz)]s,P,D \h ® h\ B (pel)) 
= <[vr + 7r _ ] Si p ijD |ji|5 _ } (vr~ |j 2 |0) or (vr~| ji\B~) ([7r + 7r~] 5i p \j 2 \0) , (6) 



since we neglect B~ annihilation contributions which are expected to be 
small [6]. Furthermore, as for the hard scattering corrections, their evalu- 
ation [9J introduces two phenomenological parameters. In Eq. ([6]) j% and 
j 2 denote the appropriate quark currents entering in Eq. Note that, 
in our approach, in the evaluation of the long distance matrix element 
([Tr + iT~]s,p,D\ji\B~) , we make the hypothesis that the transitions of B~ 
to the [tv + 'K~]s,p,d states go first through intermediate meson resonances 
Rs,p,d which then decay into a tt + it~ pair. We describe these decays by 
a vertex function modeled by assuming them to be proportional to the 
pion scalar or vector form factors or to a relativistic Breit-Wigner formula, 
respectively. For the short distance part of the decay amplitudes propor- 
tional to a combination of the effective coefficients a^{M\M 2 ) it can be 
seen that for terms coming from the first line of the right hand side of 
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Eq. (|6j) Mi = [tt + tt~]s,p,d, M 2 = tt~ and for those from the second line 
Mi = 7r~ while M2 = [vr + 7r _ ]5 i p, the [7r + 7r~]£) transition to the vacuum 
being zero with the involved bilinear quark current j 2 in Eq. . In the 
following, when M 2 = [7r + 7r~]<j p, we assume that the NLO corrections 
Vi{M<i) and P,f(M 2 ) are evaluated at the meson resonances Rs,p position. 
Here we take Rp = p(770)° and Rs = /o(980). A similar approximation 
has been applied in Refs. S] for the [Ktt]s,p states with R P = if* (892) 
and R S = K* (1430). 

Introducing the following short distance terms, with L = S,P,D and 



with i? D = / 2 (1270), 
u(i? L vr-) = \ u { ai {R L 7r-) + a u A {R L TX-) + < (^~) " [^6(^0+ 

4(R L n-)]rZ}, (7) 

^(vr-i? 5 ) = K [-2a%(ir-R s ) + o^tt" fl s )] + A c [-2ag(7r-# 5 ) 

+a^(7T-i? 5 )] , (8) 



w(vr i?p) = A u <^ a 2 (vr fl P ) - a^(vr fl P ) + - [a 7 (vr R P ) + a 9 (vr i? P )] 



one obtains, from Eqs. ©, (H|) and ([6|), the following 5"-, P- and D-wave 
matrix elements 



K (k~(pi)[k + (P2)^(P3)]s\T p \B-) = X s u(Rstt-) + Y s v(ir-R s ), 



X p (^(Pi)[K + (P2)Tr-(P3)]p\T p \B-) = X P u(Rpir-) + Y P w(tt-R p ), 



and 




(10) 



(11) 
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A p <7r-(pi)[7r + (p2)7r-(p3)]i3|r p |S-> = X D U (R d tt-). (12) 



p=u,c 



In Eq. the chiral factor r£ is given by r£ = 2m 2 K /[{mh + m u )(m u + m^)], 
and being the u and d quark masses, respectively. The long distance 
functions Xs,p,d and Y$.p, evaluated in Appendix [A] read 



Xs = ([7r + {p 2 )iT-{p3)}s\(ub) V - A \B-) (<K-{ Pl )\(du)v-A\0) 

= -yj\xs U (Ml - s 23 ) F BR s(ml) JTM, (13) 

Y S = (7T-(pi)|(J6)« J _p S |5-) ([TT + {p 2 )TT-{ P3 )] S \(dd) sc+ps \0) 

= v/f^^^^^rrM, (14) 

X P = ([ir + (p 2 )7r-( P3 )]p\(ub)v-A\B-) (ir-(pi)\(du) V -A\0) 

= N P ^- ( Sl3 - s l2 ) A BRp (ml) FT (823), (15) 

JRp 

Y P = (ir-( Pl )\(db) V - A \B-) ([tt + (p 2 )tt-( P3 )]p\(uu) V -a\0) 

= (si 3 -s 12 )F^(s 23 )Fr( S23 ), (16) 
X D = ([TT + (p 2 )7r-(p 3 )] D \(ub) v ^ A \B-) (tt-(pi)|(Ju)v-a|0> 

= F BRn {m 2 ) \P G h D(si 2 ,s 23 ) 

^2 71 V 3 m\ D - s 2 3 - im RD T(s 2 3) ' 



(17) 



The different quantities entering the above equations are discussed below. 

The 5"-wave strength parameter xs [Eq. {13])] will be fitted together 
with the correction P-wave parameter Np [Eq. (|15j . The deviation of Np 
from 1 corresponds to the possible variation of the strength of this P-wave 
amplitude proportional to /„■/ fp P [compare Eqs. (|A.7|) and (|A.19j) ] . 

Three scalar-isoscalar /q resonances, viz. /q(600), /o(980) and /o(1400), 
are present in the tttt effective mass range, m n7T , considered here. Since some 
of them are wide, like /o(600), one could have a possible Rs dependence 
in xs- The transition form factor from B to Rg, F B Rs (m%), could also 
depend on m nn . However, one expects these dependences to be weaker 
than the effective mass dependence of the pion scalar form factor, T™*(s 23 ), 
in which all these resonances are incorporated. Therefore we assume that 
Xs and F s (m%) are constant. This hypothesis will be assessed by the 
quality of the fit obtained with our model. We shall take Rs = /o(980) for 
the evaluation of F B s {mfy and we use F s (m%) = 0.13 [TO] , 
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For the pion decay constant we take f n = 0.1304 GeV [18] . The Rp 
decay constant is denoted by fp p and the 5-meson mass by Mp. Since the 
7r + 7r~ P-wave is largely dominated by the p(770) meson we choose fp p = 
f p = 0.209 GeV [9]. The quantity B = -2 (0\qq\0)/f% is proportional to 
the quark condensate. We calculate it as Bq ~ m 2 /(m n + m^). At the 
renormalization scale [i = va\,j2 we use = 4.9 GeV and m u = = 
0.005 GeV. For the transition form factor between the B meson and Rp 
state we set A BRp (m 2 ) = 0.37 [20]. 

For the Bit scalar and vector transition form factors F Brr (s) and Ff^(s), 
we use the following light-cone sum rule parametrization developed in Ap- 
pendix A of Ref. [21] , viz. 



*P(s) = 7^ (18) 

1 - S/Sq 

B _ 0-744 0-486 
Fl (S) " l-s/M*. " W^' (19) 

with s = 33.81 GeV 2 , Mp* = 5.32 GeV and si = 40.73 GeV 2 . The 
pion non-strange scalar and vector form factors T"*(s) and F^^s) will be 
discussed in the next section. Note that [22] 

rT(s) = ^-([7r + n-}s\nn\0), (20) 

with nn = — = (uu + dd) . 

The transition form factor between the B meson and the Rp state 
F BR ° (m 2 ) is not well known [23], so it will be taken as a free parame- 
ter to be fitted. The expressions of the tensor angular distribution factor 
D(s\ 2 ,s 23 ) and of the Rp mass dependence width r(s23), similar to those 
used for the f 2 (\21{)) contribution in the BABAR Collaboration Dalitz plot 
analysis [5], are displayed in Sec. IA.5l of the Appendix [Aj The expression of 
the /2(1270) coupling to tttt, Gj 2 , is also given there. 

In summary, from the S-, P- and D-wave matrix elements ()10p . (jlll and 
(|12p . we obtain the total symmetrized amplitude for the B~ — > 7r + 7r~7r - 
decay as 

Mj ym (s 12 ,s 23 ) = -j= [Mg(s 12 ) + Mg(s 23 ) + M P (s 12 )(s 13 - s 23 ) 
+ M P (s 23 )(s 13 - s 12 ) + M D (s 12 )D(s 23 ,s 12 )+M^{s 23 )D(s 12 ,s 23 )] ,(21) 
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with 



Gf 



Mp(sij) 



and 



V2 



MgM = ^= [-XsU [Ml - Sij ) F BRs (ml)u(R s n 



Ml -ml „_ 

+B B_ 7T z^Btt 



m b - m d 



F^( Sij )v(7r-Rs) 



N P ^A* Rp (ml)u(Rpir-) + if* ( Sij )w(ir- R P ) 

JRp 



(22) 



(23) 



M D {sij) 



^u{R D .-)^=F BR °( m l) 



G h 



m R D ~ s ij - imR D T(sij) 



• (24) 



For the fully symmetrized B + -4 tt + tt ir + decay amplitude we have 



M 



t ym (si2, S 23 ) = -J= [M% (Sl 2 ) + X J(s 23 ) + Atp(si 2 )(si3 - S 23 ) 



+ M P (s 23 )(s 13 - s 12 ) + Mi(s 12 )D(s23, s 12 ) + Mi(s 23 )D(s l2 , s 23 )] , (25) 
with 



Mt PD {sij) = M S PD { Sij , X u -> A;, A c -> A:, 5- -4 B+) . (26) 



3. Scalar and vector form factors 

As shown in Ref. [23] the full knowledge of strong interaction meson- 
meson form factors is available if the meson-meson interaction is known at 
all energies. The calculation of the S- and P-wave amplitudes (|22p and 
(|23p requires the values of the scalar and vector Bit, B(tttt) and pion form 
factors. The knowledge of the B — >■ tt and i? — >• [7r7r]s 5 p transition form 
factors is needed far below the Btt and B[Tnr]s,p scattering region. One 
has then to rely on theoretical models constrained by experiment, as we do 
here for the B[tttt]s form factor, using the value (see above in the previous 
section) determined in Ref. [19]. One could also use covariant light-front 
model, like that of Ref. [25J or, if available, semi-leptonic decay analysis 
results. For the Btt form factors we take the QCD light-cone sum rule 
results of Ref. [ST] recalled above in Eqs. (fTSj) and (fTUj) . The special case of 
the pion form factors is developed below. 
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3.1. The pion scalar form factor 

In the 7T7T case, the low-energy S wave being known and modeling the 
high-energy part one can rely on the Muskhelishvili-Omnes equations |26] to 
build up the pion scalar form factors. Their evaluation from these equations 
has been discussed in Ref. [27] and followed and developed in Ref. [28]. How- 
ever here, we shall use another approach, initiated in Ref. [22J and applied, 
using a different irn scattering matrix, in Ref. [lj. Extending this last work 
by introducing three channels and keeping the off-shell contributions, the 
pion scalar form factor T"*(s) entering in the S- wave amplitude Eq. (|22p is 
modeled according to the following relativistic three coupled-channel equa- 
tions 

3 

TT{s) = I%{E) + Y. R ^ E ) H ^ E )i i = l>2,3, (27) 

3=1 

with 



- J ^E,^) —^—- J±£ (28) 

where E represents the total energy, i.e., in the irn center of mass, E = 
and p is the off-shell momentum. In Eqs (|27p and (|28p . the indices i,j = 
1, 2, 3 refer to the tttt, KK and effective (2tt)(2tt) channels, respectively. The 

center of mass momenta are kj = y/s/4 — m 2 , with nil = m^, ni2 = tuk 

and 7713 = ?7i(27r)- The T matrix is the corresponding three-channel two- 
body scattering matrix. Here we use the solution A of the three-coupled 
channel model of Refs. |29l I3U| . where the effective Tn^ 2n )= 700 MeV. The 
functions Rf(E) are the production functions responsible for the formation 
of the meson pairs before their scattering. From Eqs. (|27p and (|28p one can 
check that 

3 k r 

Im Vr(s) =-J2 -j^ T U E > k J ,k l )T]*(s)9(V~s - 2m,). (29) 

This is the same unitary relation as that of the corresponding Muskhelishvili- 
Omnes pion scalar form factors constructed in Ref. [28J [see Eq. (28) therein]. 

In Eq. ([28]) the regulator function (k 2 - + K 2 )/(p 2 + k 2 ), which reduces 
to 1 on-shell (kj = p), ensures the convergence of the integral. The range 
parameter k will be fitted to data. The choice of a separable form for 
the interaction yields analytic expressions for the T matrix elements. One 
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introduces a rank-2 separable potential in the kit channel and a rank-1 
separable potential in the KK and in the (27r)(2-7r) ones. According to the 
formalism developed in Ref. [31] and applied in Ref. [29] one has for the T 
matrix elements: 



Tu(E,p,ki) = g (k 1 )t 00 (E)g (p) + g 1 (k 1 )t 11 (E)g 1 (p) + g (k 1 )t w (E)g 1 (p) 

+ 9i(h)toi{E)g (p), 

T 21 (E,p,h) = g (ki)t 2(E)g 2 (p) + gi(h)t 12 (E)g 2 (p), 

T 31 (E,p,h) = go{kx)t 03 {E)g 3 {p) + gi (kx)t 13 {E)g 3 (jp), (30) 



where 



go (h) 
9j{h) 



Air 



1 



4^F 1 

m kf + J) ' 



j = 1,2,3. 



(31) 



The parameters f3j, j = 0, 1,2,3, of the separable form of the scattering T 
matrix are given in Table 1 of Ref. [29] (fit ^4) . 

One can extend the expressions of the reduced symmetric t(E) matrix 
elements given in terms of the separable potential parameters in Appendix 
A of Ref. |31] to the case of Ref. |29] which we use here. The Yamaguchi 
form [32] of the go (p) and gi (p) (f3T|) in the T matrix elements ([30]) leads the 
following analytic expression for Tf*(s) in Eq. (I27p 



rf (a) =RUE) + RUE){[too(E)g (k 1 ) + t 01 (E)g 1 (k 1 )]g (k 1 )F 10 (k 1 )+ 
[h^gxih) + t 10 (E)g (k 1 )}g 1 (k 1 )F ll (k 1 )}+ 
R%(E) [g (fei)t 02 (E) + 9l (h)t 12 (E)]g 2 (k 2 )F 22 (k 2 )+ 
R n 3 (E) [g (fci)to3 (E) + gi (h)t 13 (E)]g 3 (k 3 )F 33 (k 3 ), 

(32) 



where 
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F w (h) 

F u {h) 
F 2 2(k 2 ) 



h,o(ki 



go(k 1 )ho{k 1 ) , 

gi{ki)hi(kiY 

h,2{k 2 ) 
52(^2)^2(^2)' 



F 33 fc 3 = 71 \h 71 \ » 33 



with 



hi(k{) — \ — ■ 9' * — 1^2,3, 
V rrii kf + «r 

/*,(**) = /*(**), (34) 



and 



d^p 1 

^jC**) = / 7^3 9j(p) / hi(p), ( 35 ) 

E - 2J P 2 + mf + *e 



where E 1 = 2y kf + mj, i = 1,2,3. The analytical expression for these 

integrals can be found in Appendix A of Ref . [31] . 

As in Ref. [22J one constraints the Tf*(s) to satisfy the low energy behav- 
ior given by next-to-leading order one loop calculation in chiral perturbation 
theory (ChPT). One writes the expansion at low s as 

r?( 8 ) + i = 1,2,3, (36) 

with real coefficients, Tf(s) being real below the tttt threshold. Using the 
expressions obtained in NLO in ChPT for the r™*(s) given in Refs. [2"2"| 153] 
one gets, 
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fii 
J l 



and 



1 fim 2 
1 + (2U 8 



+ 



P 
+ 



LI) + 8 2m ^ + 3m ' (2^6 " LI) 



P 



1 



36vr 2 / 2 16vr 2 / 2 6 i/ 2 96vr 2 / 2 V 3 i; ' ° - 



^ (^4 + ^5) 
1 



1 



16vr 2 / 2 



1 + 



64vr 2 / 2 



1 + log 



K 



mz 



192tt 2 P 



9m 2 



(37) 



dS = 



l 

v 7 ^ 



1 + 



'"7 



+8 



48tt 2 / 2 
§m 2 K + m 2 



, m; 16mi 



m - LI) 



p 



(2Li 



m 



K 



Tfi 

1 + log ~4 



in 

J2 — 



1 

j2 



~ 2 (2LI + Lg) 



64vr 2 / 2 



m 



K 



432vr 2 / 2 m 2 



64vr 2 / 2 



1 + log 



m 



K 



64vr 2 / 2 



1 + log 



mt 



(38) 



v being the scale of dimensional regularization and / = . Further- 

more for the ChPT low-energy constants, LV, k = 4,5, 6, 8, we use the recent 



determinations of lattice QCD at v = 1 GeV as given in Table X of Ref. |34j . 
For / = 92.4 MeV, we obtain d n = 1.1957, /f = 3.1329 GeV~ 2 
and / 2 ™ = 1.6719 GeV" 2 . As in Ref. [28] we assume T%(0) 
to = and we also assume fn = 0. 



, = 0.7193 
which leads 



The real production functions are parametrized as 



RUE) 



+ t?E + uj?E 2 



l + cE A 



1,2,3, 



(39) 



the fitted parameter c controling the high energy behavior. The other pa- 
rameters, af, r™ and uf are calculated by requiring that T"(s) in Eq. ([27]) 
has the low energy expansion Eq. (|36p . These nine parameters satisfy a 
linear system of nine equations displayed in Appendix [B] Their numerical 
values, depending on the value of the range parameter k [see Eq. ([25]) ]. will 
be given in Sec. [SJ 
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3.2. The pion vector form factor 

As for the scalar case one could use the Muskhelishvili-Omnes equations 
to built up the pion vector form factor. This was done in Ref. [3] for the Ktt 
vector form factor. Here, noting that the knowledge of this form factor is 
required to describe the r~ — > 7r _ 7r°i/ T decay, we shall use the phenomeno- 
logical model of the Belle Collaboration [35J. Fitting their high statistics 
data, they built the pion vector form factor i ?1 1 r ' r (s23) by including the con- 
tribution of the three vector resonances p(770), p(1450) and p(170Q). Here 
we use the parameters given in the third column of Table VII of Ref. [35J . 

4. Observables and data fitting 
4-1. Physical observables 

The symmetrized B~ — > tt^tt^^ amplitude (|2ip depends on the two 
effective tttt masses, myi = \fs\i and 77123 = of the Dalitz plot. In the 
center of mass of tt~(pi) and ir + (p2), the pion momenta fulfill the equations 



\pt\ 



12 



Am 2 , 



\fi\ 



2m\2 



Ml - (mi2 + m n ) 



M B ~ ( m i2 - m TtY 



(40) 



and the cosine of the helicity angle 9 between the direction of P2 and that 
of P3 reads 



cos 6 



1 



(41) 



For fixed values of the effective mass m\2, the variables cos# and 77123 are 
equivalent. 

The double differential B~ — > it~it + it~ branching fraction is 



d 2 B- 



dmi2 dcos9 T B 8(2vr) 3 M| 



1 m 12 \p2\\P3\ I . A - , x |2 

M sym {s 1 2,S23)\ , 



(42) 



where is the total width of the B~ . Since the Dalitz plot is symmetric 
under the interchange of 77112 and 77123, one can limit the integration range 
on 77123 t° the values larger than 77712; hence, the differential effective mass 
distribution reads 



dB~ 



dm 



12 



COS 9 a 



d 2 B- 



dm\2 dcosO 



dcos 9, 



(43) 



0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 

E(GeV) 



Fig. 1. Modulus of the pion scalar form factor T™ (solid line), obtained in our fit 
using the NLO aP i with k — 2 GeV and for which the fitted parameter c = (19.5 ± 
4.2) GeV -4 , compared to that calculated in Ref. [37] using the Muskhelishvili- 
Omncs equations (double-dash dot line). The dash-dot line (for c = 15.3 GeV -4 ) 
and the dashed one (for c = 23.7 GeV -4 ) represent the variation of the T™ modulus 
when c varies within its error band. 

where cosO g corresponds to the value of cos 6 in Eq. (|4ip with m-12 = rri23, 
viz., 

cos 6 g = {Ml - 3m\ 2 + 3m\ ) . (44) 

The variable mi2 in Eq. (|43p is also called the light (or minimal) effective 
mass Tn m i n while 717-23 is th e heavy (or maximal) effective mass, m max . The 
U - —J- 7r - 7r + 7r - branching fraction is then twice the integral of the differ- 
ential branching fraction (|43p over myi- 

4.2. Data fitting 

We aim at describing the experimental ir + ir~ distributions obtained by 
the BABAR Collaboration in the Dalitz plot analysis of the — > 7r ± 7r =l= 7r =F 
decays [5]. Two different background distributions, related to the qq and 
the BB components, are subtracted from Fig. 4 of Ref. [5]. Six light ef- 
fective 7r + 7r - mass distributions are extracted for B + and B~ decays with 
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a subdivision of the data into positive and negative values of the cosine of 
the helicity angle 9. For the B + and B~ distributions we reject two data 
points corresponding to the 7r + 7r - effective masses equal to 485 and 515 
MeV. Also two points at 470 and 530 MeV for the four mass distributions 
with cos (9 > or with cos# < are not taken into account. This is done to 
exclude the possible contribution of the decay processes 7T . 

As a by-product of the background subtraction, five data points, with a 
small number of events, have negative values with small statistical errors. 
For these five data points we increase their errors to values corresponding 
to those of the points lying in a close vicinity. This is done at 1385 MeV for 
the B~ distribution, at 1475 MeV for the B+ one, at 290 and 1610 MeV for 
the B~ distribution with cos 9 > and at 1490 MeV for the B~ one with 
cos 9 < 0. 

We perform a \ 2 fit to the 170 data points corresponding to the six 
invariant mass distributions described above. In addition, we include the 
experimental branching ratio for the B^ — > p(770)°7r =l= , p(770)° — > 7r + 7r~ 
decay channel. The theoretical distributions are normalized to the number 
of experimental events in the analyzed range from 290 up to 1640 MeV. In 
the fits, done for a fixed value of the range parameter k entering Eqs. (|28p 
[see Sec. 5], the following four parameters were varied: the production 
functions Rf(E) [Eq. (|39|) ] parameter c, the real S-wave strength parameter 
XSi the real P-wave correction parameter Np [Eq. f)15[) ] and the transition 
form factor F BRD (ml) [Eq. (fTTll ]. 

5. Results and discussion 

In the fits to the selected BABAR data as described in the previous 
section, the CKM matrix elements [see Eq. ([2])] are calculated with A = 
0.2257, A = 0.814, p = 0.135 and fj = 0.349 p2] which leads to X u = 
1.26 x 10~ 3 - i 3.27 x 10~ 3 and A c = -9.35 x 10" 3 - i 1.72 x 10~ 6 . The 
LO contributions of the Wilson coefficients to the Eq. are given in 
the second and fourth columns of Table [TJ The sum of the leading order 
coefficient plus the next-to-leading order vertex and penguin corrections for 
the coefficients, entering into n(i?s j p7r~) [Eq. Q], v(ir~ Rs) [Eq. (JSj) ] and 
w(ir~Rp) [Eq. ©], are displayed in columns three and five, respectively. It 
can be seen that the NLO corrections are relatively small except for the 
coefficient 02 which, however, has only a small contribution to the decay 
amplitude. The corrections are calculated according to Refs. [7] and [1] 
using the Gegenbauer moments for pions taken from the Table 2 of Ref. [71 
and the corresponding moments for the p meson from Table 1 of Ref. [36] . In 
the calculation of the coefficients a p G {i:~Rs) and ag(-K~Rs), contributing to 
v(ir~Rs), we apply the method explained in Appendix A of Ref. |llj . Here 



18 




Fig. 2. The ir + ir~ light effective mass distributions from the fit to the BABAR 
experimental data [5], a) for the B~ decays and b) for the B + decays. The long- 
dash line represents the S- wave contribution of our model, the dot line that of the 
P wave, the short-dash line that of the D wave and the dot-dash line that of the 
interference term. The solid line corresponds to the sum of these contributions. 
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Fig. 3. As in Fig. [2] but for the B decays a) with cos 6* < and b) with cos 9 > 0. 



the renormalization scale fj, = and we take for the strong coupling 

constant a s (mi ) /2) = 0.303. 
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Fig. 4. As in Fig. [3] but for the B + decays. 



Table 1. Leading order (LO) and next-to-leading order (NLO) coefficients 
O'iiRs.p^^), o^(n~Rs) (in parentheses) and a^{-n~Rp) [see Eq. (0] entering into 
u(Rs,pTr~) [Eq. ©], v(w- R s )\Eq. ©] and w(ir- R P ) [Eq. ©], respectively. The 
NLO coefficients are the sum of the LO coefficients plus next-to-leading order ver- 
tex and penguin corrections. Here the renormalization scale is /i = mt,/2. The 
superscript p is omitted for i = 1, 2, 3, 5, 7 and 9, the penguin corrections being 
zero for these cases. 

a%{Rs,pn~) a^(TT~R s ,p) 





LO 


NLO 


LO 


NLO 


"l 


1.039 


1.071 + i0.03 






a-2 






0.084 


-0.041 - i0.114 


a\ 


-0.044 


-0.032 - z'0.019 


-0.044 


-0.032 - i0.019 


a\ 


-0.044 


-0.039 - i0.007 


-0.044 


-0.039 - i0.007 


u 
"6 


-0.062 


-0.057 - z'0.017 


(-0.062) 


(-0.075 - t0.017) 


a% 


-0.062 


-0.062 - i0.004 


(-0.062) 


(-0.079 - i0.004) 


a 7 






0.0001 


0.0 + i0.0001 


ag 


0.0007 


0.0008 + iO.O 


(0.0007) 


(0.0007 + »0.0) 


«8 


0.0007 


0.0008 + iO.O 


(0.0007) 


(0.0007 + iO.O) 


«<) 






-0.0094 


-0.0097 - i0.0003 


"10 


-0.0009 


0.0006 + i0.0010 


-0.0009 


0.0006 + i0.0010 


a 10 


-0.0009 


0.0006 + iO.0010 


-0.0009 


0.0006 + iO.0010 



There are five free parameters at our disposal. Two of them, the reg- 
ulator range k and the high energy cut-off c of the production functions 
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Table 2. Parameters of the production functions R™(E) Eq. (IMl) for k — 2 GeV 
i g? (GeV 1 ) off (GeV- 2 ) 

1 0.7095 -0.2707 1.6251 

2 0.5759 -0.0032 1.4171 

3 1.003 0.3724 2.7427 



[Eq. (|39p ] are linked to the determination of the S-wave form factor. 
The other three, xs> Np and F BR ° (m 2 ) are related to the strength of the 
S, P and D amplitudes, respectively. The range k should be larger than 0.8 
GeV which is the on-shell pion momentum approximately equal to the half 
of the effective m 7Tn upper limit ~ 1.64 GeV which we used. In our fits we 
find that the total x 2 decreases slowly when k decreases from the high value 
of 5 GeV. Here we fix the range parameter k to be 2 GeV. We perform two 
fits for the full S + P + D-wave amplitude calculated with the NLO and with 
the LO a? coefficients. Hereafter the quoted results given inside parentheses 
correspond to the numbers obtained in the second fit. The quoted errors on 
our results come from the statistical errors in the experimental data. 

A good overall agreement with BABAR's data is achieved with c = 
19.5 ±4.2 (18.9 ±4.1) GeV" 4 , xs = -19.4 ± 2.5 (-19.8 ± 2.6) GeV" 1 , 
N P = 1.122 ± 0.034 (1.015 ± 0.035) and F BR °(ml) = 0.0977 ± 0.0070 
(0.1010 ± 0.0072). The total x 2 is equal to 231.6 (233.5) for the 171 ex- 
perimental points of the fit. For both fits the branching fraction for the 
B ± -> ^(770)°^, p(770)° -> vr+vr" decay is (8.1 ± 0.5) x 10~ 6 , to be 
compared with the BABAR Collaboration determination of (8.1 ± 0.5 ± 
1.2l?; 4 ) x 10~ 6 w (8.1 ± 1.6) x 10~ 6 from their isobar model analysis [5j. 
Note that for the LO fit we explain essentially the BABAR Collaboration's 
result without significant modification of the P wave normalization, the pa- 
rameter N P w 1.02 being close to 1. For the NLO fit, N P « 1.12 ± 0.03 
and one can compare Np — 1 ~ 25% with the average 20% error of the 
experimental branching ratio. 

The CP average total branching fraction of the — > 7r ± 7r T 7r ± decays 
calculated in the NLO fit is equal to (15.2 ±1.1) x 10~ 6 to be compared to 
the measured value of (15.2 ± 0.6 ± 1.2l°; 4 ) x 10~ 6 (table III of Ref. [3]). 
The branching fraction for the S wave equals to (2.3 ± 0.4) x 10~ 6 and that 
for the D wave is (2.8 ± 0.4) x 10~ 6 . The latter value is larger than the 
branching fraction for the / 2 (1270)vr ± , (0.9 ± 0.2 ± 0.1 j^i) x 10~ 6 , deter- 
mined in Ref. [5]. In the experimental analysis the two resonances, namely 
/2(1270) and /o(1370), overlap to a large extent, which makes their sep- 
aration difficult and some part of the branching fraction obtained for one 
resonance could have been attributed to the other one. The isobar model 
analysis of Ref. |Sj gives (2.9 ± 0.5 ± 0.5±S:b) x 10_6 for the branching frac- 
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tion of /o( 1370)71^ . Then, the sum of the branching fractions for the two 
resonances equals to 3.8 x 10~ 6 . This value compares well with the branch- 
ing fraction of 3.6 x 10 -6 obtained by integrating our distribution in the 
m n7T range between 1.0 and 1.64 GeV in which both /2(1270) and /o(1370) 
give their dominant contributions. In our model the D-wave contribution is 
dominant in this range. Let us note that the value we obtain for the transi- 
tion form factor F BRd (m^) is 29% larger than the value 0.076 given in Table 
1 of Ref. [23j for the ISGW2 model. The 5-wave contribution represents 
here as much as 15% of the total branching fraction. This contribution is 
of the same order as that of the p(1450) and p(1700) which also represents 
15 % of the total P-wave contribution. 

Before comparing our effective mass distributions to the experimental 
ones, we now give our result for the pion scalar form factor r™(s). With the 
fixed value of k = 2 GeV used in the fits, one obtains for the a", r™ and 
uif, i = 1,2,3, entering into Eq. (|39p . the values given in Tabled Then, in 
Fig. [H we show the modulus of the pion scalar form factor obtained using 
the NLO coefficients aP i for the fitted value of the parameter c = (19.5 ± 4.2 
GeV) -4 together with its envelope when c varies within its error band. It is 
also compared to that of the scalar form factor calculated by Moussallam |37| 
solving the Muskhelishvili-Omnes equations [26] with a high-energy ansatz 
starting at 2 GeV and the same low-energy three coupled-channel scattering 
T-matrix as in our model (see Sec. 3.1). However, in his calculation the 
off-diagonal matrix elements T\^{E,ki,p) and T23(ki, E,p) are set to zero 
in the unphysical region E < 2m3 = 1.4 GeV. Let us remind here that 
the imaginary parts of these two pion form factors satisfy exactly the same 
relation given by Eq. (|29p . The functional dependence of both T"*(s) moduli 
is quite similar. It can be seen in Fig. Q] that, within our model, the needed 
r"*(s) is relatively well constrained. If we fix k = 3 GeV then the fit to 
BABAR data gives c = (30.4 ±6.6) GeV" 4 , xs = (-20.2 ±2.9) GeV" 1 with 
a total x 2 °f 234.1. In this range of variation of the strongly correlated k and 
c parameters, we have checked that the scalar form factor varies smoothly. 
The corresponding values of the strength parameter xs, being very close to 
—20 GeV -4 , are not sensitive to these variations. For k = 3 GeV the values 
of the branching fractions for the different tttt waves stay within the error 
bands of those for the k = 2 GeV case. 

The threshold behavior of our pion form factor is governed by the chiral 
perturbation expansion Eq. (|36p . These ChPT constraints, not explicitly 
included in Moussallam's case, lead to r™*(s) moduli of both approaches to 
differ only slightly near the tttt threshold. Above the tttt threshold, there is 
a maximum corresponding to the /o(600) resonance, then close to 1 GeV 
a characteristic dip due to the /o(980) and finally, below the spike at 1.4 
GeV related to the opening of the third channel, there is some enhance- 



22 



ment generated by the /o(1400) present in the tttt three-channel model used 
here [29\ 130] . The third threshold energy equal to 1.4 GeV is a parameter 
representing twice the mass of the effective two-pion mass m( 27r ) used to 
account for the four pion decays of scalar mesons (see Ref. |29|). Thus, in 
nature there is no such sharp energy behavior. These characteristic features 
of the pion scalar form factor T™(s) are essential to obtain a good fit of the 
experimental effective mass distributions of the to 3tt decays. 

The results of the fit on the experimental distributions, obtained using 
the NLO coefficients aP i in the i? ± — > 7r =l= 7r =F 7r ± amplitudes, are displayed in 
Figs. El [3] and 01 The p(770)-resonance contribution dominates the it + it~ 
spectrum, but that of the S-wave is non negligible. As seen, the S'-wave 
part is sizable near 500 MeV which is related to the contribution of the 
scalar resonance /o(600), not explicitely included in the BABAR Dalitz 
plot analysis [Sj. In the 1 GeV range the /o(980) resonance is not observed 
as a peak in the ir + ir~ spectrum. This fact is easily explained in our model 
since the decay amplitudes are proportional to the pion scalar form factor 
which has a dip near 1 GeV as seen in Fig. [TJ Around 1.3 GeV there is a 
maximum coming from the contribution of the /2(1270) resonance. Near 1.4 
GeV the /o(1400) scalar resonance |29} 130] gives only a tiny enhancement 
in the distributions. 

Figure [2] exhibits a small CP asymmetry, the B~ and B + effective mass 
distributions being very close. Summing the number of experimental events 
in the m 1T + 1T - range between 290 and 1640 MeV one finds 616 events for the 
B~ decay and 606 for that of the B + . This leads to a CP asymmetry of 
(0.8±4.8)% which can be compared to the values of (1.7±0.2)% for the NLO 
fit and (—0.06 ± 0.08)% for the LO fit. Taking into account the statistical 
error of 4.8% and adding to it a few percent systematic error one sees that 
both fits agree with experiment. Let us recall here the experimental value 
of the CP asymmetry A C p = (3.2 ± 4.4 ± 3.1±l%)% for the total sample of 
7r =l= 7r =F 7r = ' = events [5]. For the particular decay mode, namely for the B decay 
into p(770)°7r ± , / o(770)° — > tt + tt~, the isobar model analysis gives Acp = 
(18 ± 7± 5±u)%, while from our model we get 3.6% ± 0.2% (-0.03% ± 
0.001%). Note here that the asymmetries obtained for the fit corresponding 
to the amplitudes calculated with the real LO aP i coefficients are quite small 
as it could have been expected. 

Figures [3] and H] show a spectacular feature, namely that the interfer- 
ence term of the S, P and D waves is quite important under the p(770)° 
maximum. Here the S-P interference dominates. The sign of this inter- 
ference term depends on the sign of cos 6, so the p peak is reduced for the 
negative values of cos 9 and enhanced for the positive values. This is a clear 
indication that the ir + Tr~ effective mass distribution cannot be reproduced 
without the 5-wave contribution. If we try to fit the data without the S- 
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wave amplitude then we obtain a poor fit with x 2 = 316.3. In this case 
the effective mass distributions are not well described below 600 MeV and 
also under the p maximum. One striking feature is that the interference 
terms allow an extremely good representation of the separate cos 9 < and 
cos 6 > spectra for the B + decays (Fig. and yield for the full spectrum 
[Fig. [2b)] a x 2 /P oint of L07 - Tne fit of the separate B~ spectra (Fig. EJ) 
is less satisfactory whereas that of the full spectrum [Fig. [2h.)] is almost 
perfect with a x 2 /P orn t of 1.2. 

6. Summary and outlook 

The present paper is a continuation of our efforts [TJ [21 El H] in con- 
straining theoretically the meson-meson final state strong interactions in 
hadronic charmless three-body B decays. If the strong interaction ampli- 
tudes are sufficiently well understood then one can improve the precision of 
the weak interaction amplitudes extracted from these reactions. 

Our theoretical model for the B^ — > tt^tt^tt^ is based on the application 
of the QCD factorization [61 \7\ El [13] to quasi two-body processes in which 
only two of the three produced pions interact strongly, forming either an 
S-, P- or L>-wave state. One assumes that the third pion, being fast in 
the 5-meson decay frame, does not interact with this pair. This hypothesis 
is mainly valid in a limited range of the 7r + 7r _ effective mass, here taken 
between the tttt threshold and 1.64 GeV. 

The short-distance interaction part of the decay amplitudes describes 
the flavor changing processes b — > uud and b — > ddd. It is proportional to 
Cabibbo-Kobayashi-Maskawa matrix elements multiplied by effective coef- 
ficients calculable in the perturbative QCD formalism. This short-distance 
amplitude is multiplied by a long-distance contribution expressed in terms 
of two products. The first one is the product of the pion decay constant by 
the B — > 7T7T transition matrix element and the second one is the product of 
the pion form factor by the B — > 7r transition form factor. The parametriza- 
tion [Eqs. (|18p . ([19p ] of the scalar and vector B to n transition form factors 
follow from the light-cone sum rule study of Ref . [21] . 

The effective Wilson coefficients are calculated to next-to-leading order 
in the strong coupling constant. They include vertex and penguin correc- 
tions but neither hard-scattering ones nor annihilation contributions since 
these last two terms contain unknown phenomenological parameters related 
to amplitude divergences [9]. We find that these vertex and penguin cor- 
rections are small in comparison to the leading order term (see Table [1]) . 
However, they allow to generate some non-zero CP asymmetries. 

We then assume the B to tttt transition matrix element to be equal to 
the product of the B to intermediate meson transition form factor by the 
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decay amplitude of this meson into two pions being either in S, P or D wave. 
The next step is to suppose the latter decay amplitude to be proportional 
to the pion non-strange scalar or vector form factor depending on the wave 
studied. For the S wave the proportionality factor is given by a fitted 
parameter xs an d for the P wave it is related to the inverse of the p decay 
constant. For the limited range of the effective tttt mass, from tttt threshold 
to 1.64 GeV, the B — > tttt transition form factors are taken as constants 
given by the B -)■ / (980) [E] and by the B -> p(770) [2DJ transition form 
factors at q 2 = m 2 . The decay amplitude for the tttt D wave is described by 
a relativistic Breit-Wigner formula and the not well known B to /2(1270) 
transition form factor is fitted. We find F Bf2 (ml) = 0.098 ± 0.007. 

The pion scalar form factor is modeled by the unitary relativistic three 
coupled-channel equation ([271) using the tttt, KK and effective (27r)(27r) 
scattering T matrix of Refs. |29[ 130] . This form factor depends on two 
fitted parameters: the first one k insures the convergence of the involved 
integrals and the second one, c, controls the high-energy behavior of the 
production functions accountable for the meson pair formation. The pion 
vector form factor takes into account the contribution of the p(770), p(1450) 
and p(1700), and follows from the parametrization of the Belle Collaboration 
in their study of the semi-leptonic r~ — > tt~tt^v t decays. For the P-wave 
amplitude we introduce a fitted correction factor Np. 

We obtain a good fit to the tttt effective mass distributions of the BABAR 
Collaboration data of the B^ — > tt^tt^tt^ decays [5]. The value of the 
branching fraction for the B ± -)• ^(770)°^ decays, (8.1 ± 0.7 ± 1.2±5;f ) x 
10~ 6 , is well reproduced with the correction factor Np close to 1. This shows 
that the QCD factorization gives the right strength of the B to pn decay am- 
plitude. The tt + tt~ spectra are dominated by the p(770)° resonance but, at 
low effective mass, the S- wave contribution is sizable. Here the /o(600) reso- 
nance manifests its presence. Furthermore one observes a strong interference 
of the S and P waves in the event distributions for cos 9 > and cos < 0. 
Here the /o(980) is not directly visible as a peak, since the pion scalar form 
factor has a dip near 1 GeV. The surplus of events in the tt + tt~ effective mass 
close to 1.25 GeV is well described by the contribution of the /2(1270) res- 
onance. The branching fraction for B^ — > /2(1270)7r :l= , /2(1270) -4 tt + tt~ 
decay is found to be of (2.8 ±0.4) x 10 -6 . At 1.4 GeV, the tiny maximum of 
the 5-wave distribution comes from the scalar resonance /o(1400) |29[ 130] . 

Our model yields a unified description of the contribution of the three 
scalar resonances /o(600), /o(980) and /o(1400) in terms of one function: 
the pion non-strange scalar form factor. This reduces strongly the number 
of needed free parameters to analyze the Dalitz plot. The functional form 
of our S- wave amplitude [Eq. (J22J], proportional to r"*(s), could be used in 
Dalitz-plot analyses and the table of T"*(s) values can be sent upon request. 
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The strong interaction phases of the decay amplitudes are constrained 
by unitarity and meson-meson data. Their determination should help in 
the extraction of the weak angle phase 7 or 4> 3 equal to arg(— A*/A*). Of 
course new experimental data with better statistics would be welcome. One 
expects — > tt^tt^tt^ 1 events from the Belle Collaboration, and probably, 
in the near future, from LHCb and from the near term super B factories. 

The authors are obliged to Bachir Moussallam for providing them the 
values of his pion scalar form factor and to Gagan Bihari Mohanty 

for useful comments on the BABAR data. We are very grateful to Maria 
Rozariska, Bachir Moussallam, Eli Ben-Haim and Jose Ocariz for helpful 
discussions. This work has been supported in part by the Polish Ministry 
of Science and Higher Education (grant No N N202 248135) and by the 
IN2P3- Polish Laboratories Convention (project No 08-127). 

Appendix A 
Long-distance functions Xs,p,d and Ys,p 



Appendix A.l The function X$ from the S"-wave amplitude 
proportional to BR$ transition matrix element 

From Eq. (|13j) the function X$ reads 

X S = ([ir + (p 2 )ir-(p 3 )} s \(ub) v „ A \B-) <7r|(<Iu)y_A|0) 

= G n Rs7T+7r -(s 23 )(R s \(ub) V - A \B-) (TT-\(du) v „ A \0) , (A.l) 

where the vertex function +7r _ (S23) describes the Rs decay into a 
[7r + 7r - ]s pair. The B to Rs transition matrix element reads (see e.g. 
Eq. (B6) of Ref. [J2]) 



(Rs(P2 +P3)\uj"(l - -y 5 )b\B-(p B )) 



mi 



A 1 
-Pi 



Ff^imD+^^p^ 

TYl— 



m„ 



(A.2) 



where Fq Rs (m^) and s (m%) are the BRs scalar and vector form fac- 
tors, respectively. The pion decay constant /„- is defined as 



(vr CPi)|rfT/x(l - 7s)«|0) = iUpip. 



(A.3) 
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The product of Eqs. (jX2]l and (TOl) yields 

(Rs\{vb) V - A \B-) (7r-|(du) V -A|0) = -(Ml - S23 )/ W i^ Rs (ml 
The vertex function G^ s7T+n _ (S23), as in Ref. [2], is modeled by 

<[vr+^-] 5 |nn|0> = G^ s7r+7r _ (s 23 ) (R s \nn\0) . 
An effective scalar decay constant /K can be introduced with 

(i?s|nn|0)=m« s /^. 
From Eqs. (fO|) . ([2Djl and HjKty one obtains 



^ S7r+7r -( S23 ) = yf X5 rr 



O23) 



^2 v^^o 
3 ™fl s Srs 



(A.4) 

(A.5) 
(A.6) 

(A.7) 



with 



V^o 



(A. 



The effective scalar decay constant has a role comparable to the Rp decay 
constant as can be seen comparing Eqs. (|A.7|) and (|A.19p . The product of 
Eqs. (|A"7T)l . (jA~2ll and ([O]) gives 

*s = ~yj\xs U {Ml - s 23 ) F BR s(ml) rr(s 23 ). (A.9) 

Appendix A. 2 The function Yg from the S-wave amplitude 
proportional to Bit transition matrix element 



From Eq. (j!4[) one has 



Yg = (7T-\(db) 



sc—ps 



B 



[K + (P2)7r-(P3)]s\(dd) 



sc+ps I 



= (7r-|J6|B-)<[7r + (p 2 )7r-(p 3 )]5|^|0>. (A.10) 
From the Dirac equations satisfied by b(pB) and d{p\) one obtains 



[tt (pi)\d(pi)b(p B )\B (p B ) 



tt (pi 



7/ x7-(PB-Pl),/ \ 

rf(Pi) bypB) 

m b - m d 



B-{p B , 
(A.11) 
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The B to 7r transition matrix element (it \(db)v-A\B ), entering into the 
above expression, can be written as (see e.g. Eq. (5) of Ref. |3j) 



(Tr-( Pl )\d^(l-^)b\B-(p B )) 
.„ Ml -ml „ 

(pb+piY- 2 V 



Ff V) + M * g2 m V iW), (A-12) 



where F^ n {q 2 ) and F^{q 2 ) are the I?7r scalar and vector form factors, 
respectively and q = Pb — Pi = P2 + £>3- Using Eqs. ()A.12j) and (f20|) in 
Eq. ([A3U]) . yields 

*S = \/|b rr( S23 ) M ^~ m " F B *(s 23 ). (A.13) 
V 3 m b -m d 

Appendix A. 3 The function A"p from the P-wave amplitude 
proportional to BRp transition matrix element 

From Eq. (|15|) one has for the function Xp (see Eq. (3.1) of Ref. 

X P = ([ir + (p 2 )Tr~(p 3 )]p\(ub) V - A \B-) (Ti~\(du) v ^ A \0) 

_ G R P n+7T-( S ^) 



V2 



■ (P2 -p 3 ) <i?p|(n6)y_A|S"> <7r-|(d«)y_A|0) (A.14) 



where the -Rp decay into a [7r + 7r ]p pair is described by the vertex function 
Gj£ p7r +„~( s 23)" Here e represents the polarization vector of the P-wave 
meson i? P . The factor l/\/2 comes from the fact that Rp represents the 
p(770)°. As seen from e.g. Eq. (B6) of Ref. H2J or Eq. (24) of Ref. [6], 



(R P (p 2 + P3)\(ub) V - A \B (p B )) = -i 2m Rp f B pi A BRp (p\) 

Pi 

+other terms. (A. 15) 

The "other terms" do not give any contribution when multiplying this ma- 
trix element by that given in Eq. (|A.3p . Plugging this expression into 
Eq. (|A.14p one has a product of polarization vectors and the sum over 
the three possible polarization eigenvalues of the state Rp should be done. 
From 

£ #P)^(p) = -C^-^)> (A. 16) 

A=0,±1 P 
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one obtains 

^2 e A • (p 2 - P3)e x " • p B = -Pi ■ (pa - P3)- (A.17) 
A=0,±1 

Then 

X P = N P ^- ( Sl3 - s l2 ) A$ R "(ml) FTM. (A.18) 
JRp 

Above, as shown in Ref. [3] for the if* (892) — > (Kir)p decay case [see their 
Eq. (D9)], we have parametrized the Rptt + tt~ vertex function in terms of 
the pion vector form factor Fi n (s 2 3). One has 

G Rpn+n ~M = N P V* JT(«23), (A.19) 
mRpJRp 

fR P being the charged Rp decay constant. Above we have introduced a 
parameter Np to take into account the possible deviation of the strength of 
the P wave, here proportional to 1/ Jr p . 

Appendix A. 4 The function Yp from the P-wave amplitude 

proportional to the Bit transition matrix element 

From Eq. (Q2D 



Y P = (TC-\(db) V -A\B-) ([7T + (p 2 )7r-(p 3 )}p\(uu) V -A\0) . (A.20) 
The pion vector form factor is defined by (see e.g. Eq. (36) of Ref. [6]) 



(Rp\(iLu) V _ A \0) = ([7T + (p 2 )Tr-(p 3 )}p\u lfl (l- l5 )u\0) = - (p 2 - P3 ) ^ ifV). 

(A.21) 

The minus sign arises from the definition of the form factor Fi w (q 2 ) which 
contains a plus sign for a (dd)v-A current [similar to Eq. (|A.12| . then as 
p° = \ j \pl{uu — dd), there will be a minus sign for a (uu)v-A current. The 
product of Eqs. (|A.12|) and (|A.21|) gives 

Y P = -2 Pl - (p2-P3)F^(q 2 )Fr(q 2 ) = («i 3 " s 12 )F^ (q 2 )Fr (q 2 )- 

(A.22) 

Appendix A. 5 The function Xp> from the D-wave amplitude 
proportional to BRp> transition matrix element 

From Eq. (|17|) one has 
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Xd = ([7r + ( P 2)7T-( P3 )}D\(ub)v~A\B-) (7r-(pi)|(du)v_A|0) 

1 2 

= 7/|GW+7r-( s 23) J] e "/3P2P3 (- R D(PU)I(^)V-A|-B _ ) 
A— — 2 

(ir-CpiJKduJv-AlO), (A.23) 

with pd = j»2 + £>3- The factor of l/v2 is due to the quark content of the 
resonance Rd [the meson /2(1270)]. The Rd decay into a [7r + 7r~]£) pair is 
described by the vertex function G RD1T + n -{s 23 ). Here e ai g(A) represents the 
polarization tensor of the /2(1270) and A is its spin projection (see Ref. [38] . 
p. 147). Taking Eq. (A3) for (7r - (pi)| (du) V - A \0) and Eq. (4) of Ref. [23] 
for the transition matrix element (R D {pD)\{ub)v -a\B~) we obtain 



X D = -^G RDn+7T -(s 23 )F BR °(ml) £ 6^(A)^6^(A«. (A.24) 

A— — 2 

To be consistent with the choice of normalization of Eq. (A2), we have 
multiplied by i the right hand side of Eq. (4) in Ref. |23] , One can show 
that (see Eqs. (7.7) and (7.8) of Ref. [35], p. 73) 



D(s 12 ,s 23 )= ^(%2^^(A)^ = l(|pt||^|) 2 -{^-^) 2 , (A.25) 

A=-2 6 

pi and p2 being the momenta of the 7T~(pi) and the vr + (p2) in the rest frame 
of vr + (p2) and n~(p 3 ). One obtains , with 77123 = \/ S2 3' 



Pi ■ P2 



V S13 - S12) 



m 



23 



4m 2 ., 



IP2I 



|P3| 



2m 



23 



M| - (777,23 + 



M 2 



(777,23 -mO ,(A.26) 



which allows to express Eq. ()A.25p in terms of s\ 2 and «23- The vertex 
function entering into Eq. (|A.23p is parametrized as being proportional to 
a relativistic Breit-Wigner resonance formula, we write 
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G Rd ^-(s23) = \~— ; J L wT^v . (A. 27) 



G f: 

where (see Ref. [38], p. 147) 



G f2 = rn f J 607VT 5 h ™ , T f27Tn = 0.848 T h (A.28) 
V q h 

and the mass-dependent width T(s23) can be expressed as (see Eq. (7) of 
Ref. @), 

Here Tj 2 is the total width of the /2(1270) resonance, mj 2 its mass and qf 2 
is the pion momentum in the f% cm. system. The Blatt-Weisskopf barrier 
form factor is given by [5] 

where the meson radius parameter rgw = 4 (GeV/c) -1 . Finally one has 

Xn = -±UF^(ml) 7| ^(^a) (A . 31) 
V2 V 3 m|j D - s 2 3 - ^i? D r(s 2 3) 

Appendix B 
Linear system of equations for af, r™ and 

The linear system of nine equations satisfied by the nine production function 
parameters a™, and cj", i = 1,2,3, is 



< + J>?i^(0) = d?, 



7 = 1 



£=0 



? + £ u^(o) +r, 

i=i v 



n dH Tl {E) 



dE 



E=0 



1 dPH^E) 
^ 2 j dE 2 



E=0 



-- o, 

//'. (B.l) 
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